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'. ABSTRACT 

Context. Constraints on the origin and propagation of cosmic rays can be achieved by studying the variation in 
. the spectral index of the synchrotron emission across external galaxies. 

' Aims. To determine the variation in the nonthermal radio spectral index in the nearby spiral galaxy M33 at a 
, linear resolution of 360 pc. 

' Methods. We separate the thermal and nonthermal components of the radio continuum emission without the 
assumption of a constant nonthermal spectral index. Using the Spitzer FIR data at 70 and 160 /xm and a standard 
dust model, we deredden the Ha emission. The extinction corrected Ha emission serves as a template for the 

• thermal free-free radio emission. Subtracting from the observed 3.6 cm and 20 cm emission (Effelsberg and the 
I VLA) this free-free emission, we obtain the nonthermal maps. A constant electron temperature used to obtain 
I the thermal radio intensity seems appropriate for M 33 which, unlike the Milky Way, has a shallow metallicity 
' gradient. 

, Results. For the first time, we derive the distribution of the nonthermal spectral index across a galaxy, M33. We 

• detect strong nonthermal emission from the spiral arms and star-forming regions. Wavelet analysis shows that 
at 3.6 cm the nonthermal emission is dominated by contributions from star-forming regions, while it is smoothly 

, distributed at 20 cm. For the whole galaxy, we obtain thermal fractions of 51% and 18% at 3.6 cm and 20 cm, 
' respectively. The thermal emission is slightly stronger in the southern than in the northern half of the galaxy. We 
I find a clear radial gradient of mean extinction in the galactic plane. 
Conclusions. The nonthermal spectral index map indicates that the relativistic electrons suffer energy-loss when 
diffusing from their origin in star-forming regions towards interarm regions and the outer parts of the galaxy. We 
also conclude that the radio emission is mostly nonthermal at R > 5 kpc in M33. 
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1. Introduction 

The problem of separating the two components 
of the radio continuum emission, free- free (ther- 
mal) and synchrotron (nonthermal) emission, dates 
back to the beginnings of radioastronomy. The 
mostly applied techniques are based on either as- 
suming a co nstant nonthermal spectral index (e.g. 
iKlein et al.|[l984) or a correlation between thermal 
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radio and infra-red (IR) emission (jBroadbent et al 
1989I) . 



Although the assumption of constant nonther- 
mal spectral index may be reasonable for global 
studies (Sect. 10), it does not lead to a feasible 
thermal/nonthermal distribution in detailed stud- 
ies. Under this simplification, it is not possible to 
investigate the origin and energy loss processes of 
the electron component of cosmic rays (CRs) within 
a galaxy. Discrete synchrotron emitting sources are 
mainly identified as supernova remnants. The syn- 
chrotron emission from supernova remnants can be 
described as power law (5^ ~ i^""") with a typical 
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Table 1. Images of M33 used in this study 



Wavelength 


Resolution 


rms noise 


Telescope 


20 cm 


51" 


70/iJy/beam area 


VLA+Effelsberg ' 


3.6 cm 


84" 


230/iJy/beam area 


Effelsbergi 


160 fj,m 


40" 


10 /^Jy/arcsec^ 


Spitzer'^ 


70 /im 


18" 


10 /ijy/arcsec'^ 


Spitzer^ 


6570A (Ha) 


2" (pixel size) 


0.3 cm~®pc 


KPNO'' 



^ Tabatabae i et al. (2007a) 

^ ffinz ct al (2004) and Tabatabaei et al. f2007a'l 
^ Hoopes fc"Walterbos (2000) 



spectral index of a„ ~ 0.6. Propagating from these 
sources, electrons suffer energy losses that steepens 
the power law spectrum. In the interstellar medium, 
the typical spectral index attributed to the emis- 
sion under the leakage loss is ~ 0.9, and under syn- 
chrotron loss and inverse Compton scattering ~ 1.1 
( Biermann These electrons further diffuse to 

the interarm regions and outer parts of spiral galax- 
ies within their life time. Hence, variations of the 
nonthermal spectral index should be distinguish- 
able particularly by comparing arms with interarm 
regions and outer parts of the galaxies. 

The thermal/nonthermal separation based on 
the assumption that the radio-IR correlation is 
due to a correlation between the thermal radio 
and IR emission is not generally correct, as non- 
thermal phenomena e.g. super massive black holes 
or energetic CRs may stimulate the IR emission. 
For instance, the possibility of heating the dif- 
fuse dust (emit ting in IR) by CRs was shown b y 
Helou fc Bicavl fl993i ) and iBicav fc Heloul (|l990l ). 
Moreover, supernovae, the ultimate sources of most 
of the nonthermal emission, explode close to the 
star-forming regions in which their proginator stars 
formed. Together with the increased magnetic field 
strength in the spiral arms, this causes a correlation 
between the thermal and nonthermal sources, mak- 
ing a less direct link between the IR and thermal 
emission. Finally, the slope of the radio-IR corre- 
lation depends on th e synchrotron spectral index 
( Niklas fc Bex3ll997h . 

Templates for free-free emission could be pro- 
vided by emission of recombination lines as they 
originate from within ionized regions, like the free- 
free emission. The Ha emission, the strongest 
Balmer line, is observationally most preferred, par- 
ticularly from nearby galaxies. Both the Ha and 
free-free emission are linearly proportional to the 
number of ionizing photons produced by massive 
stars. On the other hand, the Ha emission suffers 
from extinction by dust leading to an underestimate 
of the free-free emission if no extinction correction 



is made. Note that the emission of radio recombi- 
nation lines is extinction free and hence ideal for 
tracing the free-free emission. However, these kinds 
of emission are too weak from the diffuse ionized 
gas in external galaxies to be detected using the 
present facilites. For example, our attempt to map 
radio recombintion lines emission from diffuse re- 
gions within IC 342 and M 33 using the 6.2 cm re- 
ceiver of the Effelsberg Telescope was not success- 
ful. 

The nearest Scd galaxy, M33 (NGC 598), at a 



distance of 840 kpc (1" ~ 4pc, iFreedman et al 



19911 ) has been extensively studied at radio 
and IR wa velengths. With an in clination of 
i = 56° (|Regan fc Vogel I Il994j) . its spiral 
structure is well visible. The central posi- 



tion of M33 given by Ide Vaucouleurs fc Leach 
()l98ll) is RA(1950) = 1''31™01.57" and 
DEC(1950) = 30°24'15.0". The position angle of 
the m ajor axis is PA~ 23° |Deul fc van der HulstI 
19871) . So far, extinction studies in M33 have mostly 
focused on HII regions using either the Balmer (and 



Paschen) line r atios (e.g. K witter fc AUer 198ll: 



s (e.g 
i; I 



Melnick et al. 1987t [Petersen fc Gammelgaard 



19971) or the ratio of Ha to radio emission (e.g . 
Israel fc KennicuttI Il980t iDevereux et al.l 119971) . 



In the line ratio method the emission from the 
diffuse ionized medium can hardly be detected in 
most of these recombination lines. The method 
based on the radio emission only works well when 
the thermal component of the radio continuum 
is independently known. For instance, the radio 
emission from single HII regions (and not HII 
complexes with a possible nonthermal emission 
contribution like NGC 604 and NGC 595, iDodorico 



19781 : [Gordon et al.l |1993( ) may be considered as 



the thermal emission. 

Using the ISOPHOT 60 and 170 /im data. 



Hippelein et al ] ([2003) found an anti-correlation 



between the flux density ratio of Ha / 60 /im and 
170 /im flux density, suggesting the extinction to be 
related to the cold dust in M33. 
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We obtain the distribution of the dust opti- 
cal depth at the Ha wavelength for the whole 
M33 (an extinction map) by analysis of dust emis- 
sion and absorption using the high sensitivity and 
resoluti on Multiband Imaging Photometer Spitzer 
(MIPS, iRieke et alJ 12004 ) FIR data at 70 and 
160 /xm. This leads to an Ha map corrected for ex- 
tinction (de-extincted Ha map), our free- free tem- 
plate. The thermal and nonthermal maps at 3.6 
and 20 cm are obtained at an angular resolution 
of 90" (equivalent to a linear resolution of 360 pc) 
by both the new (using a free-free template) and 
the standard method (assuming a constant nonther- 
mal spectral index) and the results are compared. 
Further, we determine variations of the nonthermal 
spectral index across M33 as detected by the new 
method and discuss the exponential scale lengths of 
both thermal and nonthermal emission. 

In Sect. 3 we derive distribution of dust color 
temperature. This temperature and the 160 fiia flux 
density are used to obtain the dust optical depth 
(extinction) and its distribution in Sects. 4 and 5. 
We correct the Ha emission for the extinction and 
then conver t it to the t hermal radio emission, fol- 
lowing JDifik^so^^t^D ([2003) (Sect. 6). The non- 
thermal intensity and spectral index maps are pro- 
duced in Sects. 7 and 8. We discuss and compare 
the results from the new and standard methods 
in Sect. 9. Finally, conclusions are presented in 
Sect. 10. 



2. Data 

The radio continuum data at 3.6 and 20cm are 

I20fl7bh 



presented in Paper I ()Tabatabaei et al.ll2007bn . At 
3.6 cm, M33 was observed with the 100-m Effelsberg 
telescope of the MPIfR0. The 20 cm data were ob- 
tained from observations with the Very Large Array 
corrected for missing spacing using the 
Effelsberg data at 20 cm. We also use the latest 
combined epoches of MIP S Spitzer data at 70 and 



160um (as presented in iTabatabaei et al.ll2007a , 



2005f) . The Ha map i s from Kitt Peak Nationa l 
Observatory (KPNO) (|Hoopes fc Walterbodlioooh . 
Table [1] summarizes the data used in this work. 



^ The 100-m telescope at Effelsberg is operated by 
the Max-Planck-Institut fiir Radioastronomie (MPIfR) 
on behalf of the Max-Planck-Gesellschaft. 

^ The VLA is a facility of the National Radio 
Astronomy Observatory. The NRAO is operated by 
Associated Universities, Inc., under contract with the 
National Science Foundation. 



3. Temperature distribution of dust 

We derive the color temperature of the dust, T, 
between 70/im and 160/ini from the formula 



F2 



Bi{T) 
'B2{T) 



(1) 



Here F denotes the measured flux, B{T) the 
Planck function, v the frequency, and 1 and 2 refer 
to the wavelength of 70/im and 160/im, respectively. 
The power law index of the absorption efficiency at 
FIR wavelengths, /?, is set to 2, which should be 
appropriate for inte rstellar grains (jAnd ricssc 1 974 
Draine fc Le'3ll984l ). 

The MIPS FIR maps at 70 and 160 /im have 
been smoothed to an angular resolution of 40" 
and normalized to the same grid and center po- 
sition. The resulting dust temperature distribution 
is shown in Fig. [1] indicating variations between 19 
and 28 K. Fig. [5] displays the relative frequency of 
occurence of the temperature intervals in 10" x 10" 
pixels. The maximum is attained at 21.5 K, close 
to the mean value of 21.6 K. Warmer dust with 
T > 25 K dominates in star- forming regions and in 
the center of the galaxy. The highest temperatures 
are found in the HII complexes NGC 604, NGC 592, 
and IC 133. 



4. Optical depth at 160jum 

Having determined T, we obtain the optical depth 
from the equation 



Iimti-m — BiQOfj.m{T) [1 



(2) 



where Iieofim is the intensity. The distribution 
of Tieo^im over the disk in M33 is plotted in 
Fig. 3. Because the temperature variations are very 
moderate, the optical depth Tigo^^m usually fol- 
lows the 160Mni intensi ty map quite closely (see 
ITabatabaei et al" 2007al) . Nevertheless, for a fixed 
intensity, a warm region (26 K) has a three times 
lower optical depth than a cold one (19 K). 



5. Distribution of extinction 

To convert rigo/jm into the dust optical depth at the 
wavelength of the Ha line, thq, we have to multi- 
ply it by KHct/KieoAim, the ratio of the dust extinc- 
tion coefficient per unit mass at the correspond- 
ing wavelengths. Using the extinction curve given 
by a standard dust model for the di ffuse medium 
(see, e.g. Figure 12.8 of Kriigel 2003f l. we estimate 
THa — 2200ri6o^un. Therefore, after multiplication 
with 2200, Fig. [3] gives also the distribution of thq 
across the galaxy at a linear resolution of 360 pc. 
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Fig. 1. Dust temperature map of M33 obtained 
from l70/jm/Ii60Aim ratio (only pixels with intensity 
above 3a level were used). The angular resolution 
of 40" is shown in the lower left-hand corner of the 
map. The bar at the top gives the dust temperature 
in Kelvin. 




22 24 26 

Dust Temperature (K) 



Fig. 2. Histogram of the dust temperature map. It 
shows the population of pixels as a function of the 
temperature intervals. The number of bins used for 
this plot is 80. 



thq is around one half in the extended central 
region and in the two main arms, IN and IS, and 
it is somewhat smaller in other arms (mostly be- 
tween 0.2 and 0.4). In the center of the galaxy and 
in massive star forming regions, specifically in the 




34 30 00 33 30 

RIGHT ASCENSION (J20001 



Fig. 3. Distribution of the dust optical depth at 
160 iim wavelength in M33. The bar at the top 
shows 10^ X Tieo/im. The main northern (IN) and 
southern (IS) spiral arms are indicated. The angu- 
lar resolution of 40" is shown in the lower left-hand 
corner of the map. 



southern arm IS and the HII complexes NGC604, 
NGC595, and B690, exceeds 0.7. The highest 
dust optical depth one finds in the center of the 
galaxy and in NGC604 where thq ^ 0.97 and 0.88, 
respectively (at the linear resolution of 360 pc). 
Therefore, M33 is generally almost transparent for 
photons with A ~ 6560A propagating towards us. 

To estimate how much the detected Ha radia- 
tion has been attenuated, we note that Ha photons 
are usually emitted from sources within the galaxy. 
The optical dept h THn- therefore only giv es an upper 
limit. Following iDickinson et all ( 20031 ) . we set the 
effective thickness to Tes = fd x Tua with < 1; 
the attenuation factor for the Ha flux is then e""^""'. 

At 360 pc resolution, one may assume that the 
Ha emitters, ionized gas in HII regions and dif- 
fuse gas, are uniforml y mixed with the dust, which 
would imply fd — 0.5. [Dickinson et al.l ( 2003 ) found 
fd — 0.33 for the Milky Way (because the z- 



distribution of the Ha emission is smaller than that 
of the dust) and we adopt their value also for M33. 
But, as will be shown in Sect. 9.4, the determina- 
tion of the thermal fraction of the radio emission is 
not very sensitive to the particular choice of fd- 
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Fig. 4. Radial distribution of the mean effective op- 
tical depth at Ha {tbs) in rings of width of 0.5 kpc 
in the galactic plane for fd = 0.33, 0.5 and 1.0. The 
errors are smaller than the size of the symbols. 



Of course, it would be preferable not to 
use a uniform value fd for the whole galaxy, 
but one that is adapted to the geometry (well 
mixed diffuse rnediuin or shell- like in HII regions. 



Witt fc Gordm] |2000[ ) and the dust column den- 
sity. This needs to specify the location of the stel- 
lar sources and the absorbing dust and to solve the 
radiative transfer problem with massive numerical 
computations. 

An interesting related question is how the ex- 
tinction in M33 changes with galactic radius R. 
To investigate this, we integrated Spitzer FIR flux 
densities at 70 and 160 /im in rings of 0.5 kpc 
width in the ga lactic plane (inclination of 56°, 
Regan fc Vogel l il994,). We first derived the mean 
dust temperature of each ring and then the mean 
dust optical depth T(.g. Fig. |4] shows Tcff versus the 
galactocentric radius R. For fd = 0.33, the depen- 
dence can be described by rH„(R) = (—0.009 ± 
0.002) R + (0.24 ± 0.03). As fd rises, the slope gets 
steeper. The case /d = 1 (for the full layer), which is 
appropriate for pure background sources of M33, is 
shown for comparison. Generally, a decrease of the 
extinction with galactocentric radius is expected as 
it reflects the decrease in surface density towards 
the periphery. 

From a comparison of H a data with radio 
flux d ensities at 6.3 and 2 1 cm, Israel fc KennicufrO 
(|1980[ ) and iBerkhuiisenI (|l983 
rather steep radial gradient of 
consid ering 8 bright HII regions 
(|1982^) derived the relation tho = AHa/1-086 = 
(-0.038 ± 0.003) R + (1.45 ± 0.05), with R in ar- 
c min. Qualitatively similar results were reported 
bv lPetersen fc GammelgaardI ( 19971). using Balmer 



1983h found a 
the extinction: 



Berkhuiisen 



and Paschen emission lines, and bv lHippelein et al 



(j2003[ ) based on a study of the Ha to the 60 /im 
ISOPHOT flux r atio in HII regions. Th e only dis- 
crepant result by Devereux et al. ( 19971 ) who, from 
their pixel-to-pixel comparison of Ha and 6 cm ther- 
mal ra dio maps, did not find a ny systematic decline 



in Ay, iHippelein et al.l (|2003[ ) were able to explain 



by assuming that the radial gradient of the cold 
dust (170 Aim ISOPHOT) is stronger in the HII re- 
gions than in the diffuse gas. 

6. Determination and distribution of 
free-free emission 

From the observed Ha intensity, /, and the effective 
extinction. Toff, we derive the intrinsic Ha intensity, 
Iq, according to 



(3) 



Integration of the Ha maps out to a radius of 

7.5 kpc yields a ratio of corrected to observed total 
Ha flux density of 1.13 for fd — 0.33 and 1.25 for 
fd = 0.5. Thus only 13% (25%) of the total Ha 
emission is obscured by dust. The small extinction 
within M33 was predicte d by the wavelet study o f 

3.6 cm and Ha emission ( Tabatabaei et al. I l2007ah . 
The emission measure ( EM) follows from the H a 

intensity via the expression IValls-GabaudI (Il998h 



iHa - 9.41 X 10-**T^ 



:i-""io- 



EM 



(4) 



where the electron temperature, Te4, is in units of 
lO'^K, EM in cm~^ pc, and it is assumed that the 
optical depth of HI resonance lines is large (usu- 
ally denoted as case B). The emission measure is 
related to the continuum optical thickness, Tc, of 
the ionized gas by 



Tr = 8.235 X lO^^aT: 



-1.35, 



^GHz (1 + 0.08) EM, (5) 



with a ~ 1 (|Dickinson et all l2003h . The factor 
(1 -I- 0.08) takes into account the contribution from 
the singly ionized He. The brightness temperature 
of the radio continuum emission, Tf,, then follows 
from 



n = Te(l - e-"=) 



(6) 



Whereas the electron temperature in the Milky 
Way is known to increas e with galactocentric 
radius (I Shaver et al.l Il983l ) as a r esult 



i Shaver et al 
creasing metallicity 



Panagia 



of de- 
M33 does 



not show signiflcant variations in metallicity (e 
Willner fc Nelson-Patell2002l : iMagrini et al.ll200 
Crockett et al. ( 2006h derived the electron temper- 
ature from forbidden line ratios in 11 HII regions 
with galactocentric distances from 1 to 7 kpc. Their 
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Fig. 5. Thermal (leff) and nonthermal (right) maps at 3.6 cm. The grey-scale gives the flux density in 
/xJy/beam. The angular resolution is 90" (shown in the lower left of the images) with a grid size of 10". 
The bright HII complexes NGC604 (RA = 34™ 32.9'' & DEC = 30° 47' 19.6") and NGC 595 (RA = 
33™ 32.4" & DEC = 30° 41' 50.0") are clearly visible in both maps, whereas IC133 (RA=l'' 33™ 15.3* 
and DEC = 30° 53' 19.7") appears stronger in the nonthermal map. 

Table 2. Thermal flux density and thermal fraction Table 3. Thermal flux density and thermal fraction 



at 3.6 cm (/d 


= 0.33 and Te 


= 10000 k). 




at 20 cm (/^ 


= 0.33 and 


= 10000 k). 




Object 


Observed 


Thermal 


Thermal 


Object 


Observed 


Thermal 


Thermal 




flux density 


flux density 


fraction 




flux density 


flux density 


fraction 




(mJy) 


(mJy) 


% 




(mJy) 


(mJy) 


% 


NGC604 


50.44 ±0.13 


38.55 ±0.08 


76.4 ±0.3 


NGC604 


62.75 ±0.21 


45.46 ±0.09 


72.4 ±0.4 


NGC595 


17.95 ±0.18 


12.13±0.11 


67.6 ±1.3 


NGC595 


20.50 ±0.22 


14.31 ±0.14 


69.8 ±1.4 


IC133 


9.08 ±0.11 


2.66 ±0.04 


29.3 ±0.8 


IC133 


5.35 ±0.07 


3.11 ±0.03 


58.1 ±0.9 


B690 


4.51±0.11 


3.11±0.10 


68.9 ±3.9 


B690 


5.80 ±0.12 


3.69 ±0.12 


63.6 ±2.4 


B61/62 


3.72 ±0.11 


2.24 ±0.07 


60.2 ±3.7 


B61/62 


4.00 ±0.21 


2.69 ±0.08 


67.2 ±4.0 


IC132 


4.89 ±0.05 


3.25 ±0.02 


66.5 ±1.0 


IC132 


5.52 ±0.05 


3.81 ±0.02 


69.0 ±0.7 


IC131 


4.59 ±0.07 


3.26 ±0.05 


71.0 ±2.2 


IC131 


4.12 ±0.07 


3.84 ±0.06 


93.2 ±2.1 


NGC588 


4.64 ±0.09 


3.24 ±0.06 


69.8 ±2.6 


NGC588 


4.67 ±0.17 


3.74 ±0.06 


80.0 ±3.2 


IC142 


3.51 ±0.08 


2.73 ±0.10 


77.8 ±4.6 


IC142 


4.14 ±0.08 


3.09 ±0.13 


74.6 ±3.4 


B691 


4.60 ±0.13 


2.35 ±0.06 


51.1 ±2.6 


B691 


7.40 ±0.20 


2.43 ±0.04 


32.8 ± 1.0 


NG(:r)92 


1.28 ±0.08 


2.r()±().()6 


61..5±2.() 


K(;Cr)92 


:).l-i±0.11 


:i.27±().()8 


()().2±2.1 


M33 (mJy): 








M33 (mJy): 








R< 7.5 kpc 


761 ± 63 


391 ± 74 


51.4±4.2 


R< 7.5 kpc 


2722 ± 60 


478 ±85 


17.6 ±0.4 



Te values range from 7300 K to 12800 K with a mean 
value of 10000 K and no clear radial gradient. As 
there are no Te measurements for the diffuse ionized 
gas in M33, we adopt a fixed value of Te = 10000 K. 



The conversion factors from brightness tem- 
perature (K) to the observed radio flux density 

(Jy/beam) are 2.6 and 74.5 at 3.6 and 20 cm, re- 
spectively. The final free-free maps are shown in 
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Fig. 6. Thermal (left) and nonthermal {right) maps at 20 cm. The grey-scale gives the flux density in 
/zJy/beam. The angular resolution is 90" (shown in the lower left of the images) with a grid size of 10". 



Figs. and [H At both wavelengths, the strongest 
thermal emission emerges from HII regions, in par- 
ticular, the HII complexes NGC604 and NGC595, 
but the southern arm IS (Paper I) and the center 
of the galaxy are also very bright. At both wave- 
lengths, the thermal fraction of the diffuse emission 
is < 25% in interarm regions. The average error in 
the thermal fraction obtained using the error prop- 
agation method is 7%. 

Integrating the thermal maps in rings around 
the galaxy center out to a radius of 7.5 kpc, we ob- 
tain the total thermal flux densities and thermal 
fractions at 3.6 and 20 cm (see Table[5]and[3]). These 
Tables also give the thermal flux densities and ther- 
mal fractions at the position of the 1 1 brightest HII 
regions (the coordinates of these HII complexes are 
listed in Paper I). 

7. Distribution of nonthermal emission 

Subtracting the maps of the thermal emission from 
the observed maps at each wavelength, maps of the 
nonthermal emission are obtained (Figs. [5] and [6]). 
The latter maps exhibit diff'use emission extending 
to large radii. They also show strong features in 
the spiral arms and central region of the galaxy. 
The strongest nonthermal emission emerges from 



the HII complexes NGC604 and NGC595 at both 
3.6 and 20 cm. Typically, HII complexes host tens of 
young 0/B stars, many of which end as supernovae 
whose remnants contribute to a mixed (thermal and 
nonthermal), flat spectrum of the total radio emis- 
sion from these regions (as discovered in NGC 604 
and NGC 595, B Qdoric_o 1971; Gordon et al.lll993l: 
lYang et all 19961 ) . Supernova remnants with central 
energy sources in the form of young pulsars (Crab- 
like remnants) also have flat nonthermal spectra. 

Another strong feature in the 20 cm nonther- 
mal map belongs to the active star-forming region 
in the central southern arm IS. This confirms our 
previous conclusion about the spatial coupling of 
the nonthermal emission with the starforming re- 
gions (Paper I). There are also other point-like 
sources which are not resolved. Looking at the 
3.6 cm map, the giant HII region IC133 (RA = l'' 
33™ 15.3" and DEC = 30° 53' 19.7") is also strong, 
but not at 20 cm. Hos ting the two stronge st opti- 
cally thick HII regions (| Johnson et al.ll200l[ ). IC133 
has an inverted spectrum. The appearance of this 
source in the 3.6 cm nonthermal map is due to 
the assumption that the free-free emission is op- 
tically thin (equivalent to the thermal spectral in- 



dex of 0.1 when S 



th 



-0.1 



). Furthermore, IC133 
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Fig. 7. The wavelet spectra of the 3.6 and 20 cm 
thermal and nonthermal emission at 90" resolu- 
tion. The data points correspond to the scales 90, 
120, 160, 214, 286, 381, 509, 679, 907". The spectra 
are shown in arbitrary units. 100" corresponds to 
400 pc in M33. 



probably contains some non thermal emission as 
Schulman fc BregmanI ( 1995h found it associated 
with a bright X-ray source located in a hole in the 
HI layer of the galaxy, indicating energetic stellar 
winds and supernovae from massive stars. 



Tabatabaei et al.l (|2007ar ) discussed the wavelet 



spectrum of the total radio continuum maps. Here, 
we present the wavelet energy of the thermal and 
nonthermal emission at different scales. Fig.[7]shows 
that the wavelet spectrum of the nonthermal emis- 
sion is smoother than that of the thermal emission 
at both wavelengths. Furthermore, the distribu- 
tion of the nonthermal wavelet energy is smoother 
at 20 cm than at 3.6 cm. Scales smaller than the 
width of the spiral arms (~ 400" or 1.6 kpc) and 
larger than the size of the giant starforming regions 
(~ 100" or 0.4 kpc) at 20 cm are not as prominent 
as those at 3.6 cm, where the synchrotron photons 
have higher energies. In fact the synchrotron emis- 
sion distribution reflects the combined distribution 
of the interstellar magnetic field and the cosmic 
ray electrons. Hence, assuming that the distribu- 
tion of the interstellar magnetic field is the same 
at both wavelengths, one can conclude: the higher 
the energy of CRs, the more dominant the localized 
structures (which can be linked to the cosmic ray 
'sources' e.g. supernova remnants). This is expected 
as energy losses of CR electrons increase with elec- 
tron energy and field strength. Studying propaga- 
tion eff ects with self-consistent Ga lactic wind simu- 
lations, Breitschwerdt et al. ( 2002f l predicted a sim- 
ilar energy-dependence for the CR nucleon distri- 
bution in our Galaxy. They expect that the high 
energy (TeV) 7-rays from CR sources dominate the 




34 30 00 33 30 

RtQHT ASCENSION (J2000) 



Fig. 8. The same 20 cm nonthermal map as shown 
in Fig. [6] with supernova remnants (crosses) super- 
imposed. The supernova remnants with a signal- 
to-noise ratio of l arger than 0.7 were s elected from 



the catalogue of iGordon et al.1 (|1998( ). The circle 



sho ws an HII region wit h nonthermal spectral in- 
dex (|Gordon et al.lll999l ). 



diffuse 7-ray emission, while the Galactic 7-ray ob- 
servations in the GeV range (with EGRET and 
COS-B) have shown a roughly uniform distribution 
of the 7-ray emissivity in the Galactic plane. 

The dominant scales of the wavelet spectra of 
the thermal (at 3.6 and 20 cm) and the 3.6 cm non- 
thermal emission are the same, indicating that the 
distribution of the cosmic ray sources is similar 
to that of the thermal sources. At scales larger 
than a = 214 ", the 3.6 cm nonthermal spectrum de- 
creases twice as fast as the 20 cm nonthermal spec- 
trum. The smaller energy losses of the CR electrons 
with lower energies (at 20 cm) than of those with 
higher energies (at 3.6 cm) may cause a smoother 
distribution of the emission at all scales. This in- 
dicates that besides the CR sources embedded in 
the star-forming regions, there is another compo- 
nent of the synchrotron emission in the form of a 
diffuse disk or a galactic halo that is better visi- 
ble at 20 cm than 3.6 cm. This could be verified by 
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Fig. 9. Nonthermal spectral index (a„) map ob- 
tained from the nonthermal radio fluxes at 3.6 and 
20 cm. The spiral arms are indicated by contours 
of the total radio emission at 3.6 cm superimposed. 
Contour levels are 1.5, 4.5, and 12mJy/beam. The 
angular resolution is 90" with a grid size of 10". 
Bottom: A histogram of the nonthermal spectral 
index map. The number of bins used for this plot 
is 50. 

multiwavelength observations of edge-on galaxies at 
resolutions higher than those of the existing data. 

Fig. [HI sho w s the the supernova remnants from 
Gordon et al.l ( 1998 ) superimposed on the non- 
thermal 20 cm map. There is general coincidence 
between nonthermal features and supernova rem- 



nants, the_£OwerfiJ sources of the relativistic elec- 
trons. Gordon et al. ( 19991 ) diagnosed 17 supernova 
remnants embedded in HII regions. They also found 
some 30 HII regions with nonthermal radio compo- 
nents of which it was not clear whether they belong 
to these regions or were external radio sources. The 
circle in Fig. [8] shows the position of the strongest 
(S20 = 7.0 ± 0.2 mJy, a ~ 0.2) HII region with non- 
thermal emission. As this HII region is not resolved 
in our thermal maps with 90" resolution, it is not 
listed amongst the bright HII regions in Tables 2 
and 3. 



8. Nonthermal spectral index 

From the nonthermal radio fluxes at 3.6 and 20 cm, 
we obtained the spectral index of the nonthermal 
emission which was only computed for pixels with 
flux densities of at least three times the rms noise a 
at both frequencies. Fig. [5] shows that the nonther- 
mal spectral index, Q!„, has a clumpy distribution. 
Note that the fainter regions have steeper spectra. 
The most probable value of Q!„ distributed across 
the galaxy is 0.95 (see the lower panel in Fig. [9]). 
In the star-forming regions, the nonthermal spec- 
trum is relatively flat with an average value of a„ 
of 0.6±0.1, the typical spectral index of supernova 
remnants, but increases to 1.2 ± 0.2 in the in- 
terarm regions and outer parts of the galaxy. This 
indicates energy losses of the relativistic electrons 
while they diffuse away from their origin in star- 
forming regions towards the interarm regions and 
the outer parts of the galaxy. For the first time, 
a nonthermal spectral index map can be used to 
achieve more realistic models for the propagation 
of CR electrons. 

Integration of the nonthermal intensity in rings 
of 0.5 kpc in the galactic plane yields the mean 
nonthermal spectral index in each ring. The radial 
variation of the ring mean spectral index is shown 
in Fig. I10( where the ring mean spectral index of 
the total radio continuum emission is also plotted 
for comparison. Up to R~ 4 kpc, there is no in- 
creasing trend for the nonthermal spectral index 
with radius as for the total spectral index. Here the 
ring mean nonthermal spectral index varies between 
0.65 and 0.9, indicating that CRs are injected by 
sour ces related to ma ssive stars, and under leakage 
loss ( Biermann 19951) . Towards the outer parts of 
the galaxy the nonthermal and total spectral indices 
converge, confirming that the total radio continuum 
emission is mostly nonthermal at R> 4.5 kpc. This 
region corresponds to the synchrotron and inverse 
Compto n. loss dominated regime as a„ ~ 1 on the 
average (|Biermannlll995i ). 



10 



Tabatabaei et al.: High-resolution radio continuum survey of M33 



O nonthermo 
A observed 



1 ,2 : 
1 ,0 : 

0.6: ,A-"^' 
0.2 : 

0.0 E : : ^ , ^ L 



-A - -« 



4 

R (kpc) 



Fig. 10. Mean nonthermal spectral index in rings 
of 0.5 kpc width in the galactic plane versus galac- 
tocentric radius. The total (observed) spectral in- 
dex is also shown for comparison. The minima at 
3 <R< 4 kpc is due to the HII complex NGC604. 



9. Discussion 

9.1. Comparison with the standard metliod 

In this section, we first obtain the distribution of the 
thermal and nonthermal emission assuming that 
the nonthermal spectral index is constant across 
the galaxy (standard method). Then we compare 
the results from the two methods. 

Because the thermal emission is weak in the 
outer parts of the galaxy, one may consider the to- 
tal spectral index from these par ts as the pure non- 
therm al spectral inde^ll, Q!„ (e.g. Berkhuiisen et al 



20031) . The total spectral index map of M33 gives 



an = 1.0 ± 0.1 (see Fig. 11, iTabatabaei et al 
For a total spectral index, a, obtained from 
the observed flux densities at frequencies and 1^2 
and the constant value of a„, the thermal fraction 
at frequency vi is given by 



'^^-o.l 



([Klein et al.lll984f ). Then the thermal flux density 

th ' 

(8) 



at frequency vi , S'^l , is obtained from 



and the nonthermal flux density at frequency i^i, 
is 



QUI 



(9) 



Using the data at 3.6 and 20 cm in the above for- 
mulae, the corresponding thermal and nonthermal 
maps are derived; those at 20 cm are shown in 

^ This is confirmed by the new method, as mentioned 
in Sect. 8. 



Fig. [m The main difference between the standard 
and our new method (Fig. (5]) concerns the distri- 
bution of the nonthermal emission, while the ther- 
mal maps show almost the same structures. The 
nonthermal emission from the standard method is 
weaker than that from the new method and hardly 
shows emission from star-forming regions in the 
arms. Often it is even weaker in the arms than in be- 
tween the arms and in the outer parts of the galaxy. 

In contrast to what is assumed in the standard 
method, large variations of the nonthermal spec- 
tral index are found across M33 by the new method 
(Fig. ^ . We interpret this as clear evidence that CR 
electrons suffer energy losses diffusing away from 
their places of origin in the arms towards interarm 
and outer regions. 

Table 4 lists the thermal fractions of the 11 
bright HII regions and of M33 (from the integrated 
flux density maps up to R = 7.5 kpc) at 3.6 cm ob- 
tained from both methods. The thermal fractions 
of all the sources are larger when obtained from 
the standard method than from the new method. 
Some even exceed 100%, causing very weak or nega- 
tive nonthermal fluxes. The thermal fraction of M33 
from the standard method is 23% ± 14% higher 
than that from the new method. In Paper I, we 
used the standard method with a„ — 1.0 ± 0.1 to 
estimate the thermal fraction from the integrated 
spectrum, based on data at 35.6, 21.1, 17.4, 11.1, 
6.3, 6.2, 3.6, and 2.8 cm, as 0.49 ±0.15 which agrees 
with the value 0.51 ±0.04 obtained from the new 
method. This indicates that the assumption of a 
constant nonthermal spectral index is reasonable to 
estimate mean values for global studies, when the 
integrated spectrum is used. 



9.2. Radial scale lengths 

Fig. HH shows the exponential distributions of the 
thermal and nonthermal intensities at 3.6 and 20 cm 
with galactocentric radius. The radial profiles of 
the nonthermal emission are smoother and flatter 
than those of the thermal emission at both wave- 
lengths. The fluctuations in the thermal proflles at 
2 < R < 4 kpc cause similar fluctuations observed 
in the proflle of the to tal emission (presented in 
ITabatabaei et al.ll2007bl ). 

The same proflles obtained from the standard 
method are shown in Fig. [13] for comparison. The 
nonthermal radial proflles are slightly smoother (es- 
pecially at 3.6 cm) than those from the new method, 
caused by the assumption of a constant nonthermal 
spectral index. Although the thermal radial proflles 
exhibit the same variations as those in Fig. [T^ the 
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Fig. 11. Thermal and nonthermal 20 cm maps obtained from the standard method. The grey-scale gives 
the flux density in /zJy/beam. The angular resolution is 90" with a grid size of 10". 
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Fig. 12. Radial profiles of the flux densities of the thermal and nonthermal emission at 3 cm (left) and 
20 cm (right) obtained from the new separation method. The thermal flux densities show a bump at 
3 <R< 4 kpc which is due to the HII complex NGC604. 



over-estimation of the thermal emission from the 
standard method is obvious. 

The exponential scale length, I, is obtained by 
fitting an exponential function of the form /(R) = 
/oea;p(— R//). where Iq is the intensity at R = 0. 
Table 5 shows the exponential scale lengths of the 
thermal emission, Ith, and nonthermal emission. In, 
from both methods. Generally, the scale lengths of 
the nonthermal emission are larger than those of the 



thermal emission (by a factor of ^ 2). In Sect. 7, we 
show that the CR sources follow the distribution 
of the thermal sources, i.e. star-forming regions. 
Hence, the radial distribution of the CR sources 
decreases faster than that of the synchrotron emis- 
sion. This is a direct observational result indicating 
diffusion of the cosmic rays from their places of ori- 
gin to larger distances. 
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Fig. 13. Radial profiles of the flux densities of the thermal and nonthermal emission at 3 cm (left) and 
20 cm (right) obtained from the standard separation method. The thermal flux densities show a bump at 
3 < R < 4 kpc which is due to the HII complex NGC604. 

Table 4. Thermal fractions at 3.6 cm obtained from the new and standard methods. 



Object 


"(3.6,20) 


Fth (new 


Fth (standard 






method) % 


method) % 


NGC604 


0.12±0.01 


76.4 ±0.3 


99.2 ±0.4 


NGC595 


0.07 ±0.03 


67.6 ±1.3 


101.6 ±1.4 


IC133 


-0.30 ±0.03 


29.3 ±0.8 


106.5 ±1.6 


B690 


0.14 ±0.05 


68.9 ±3.9 


99.1 ±3.6 


B61/62 


0.04 ±0.08 


60.2 ±3.7 


104.3 ±4.2 


IC132 


0.07 ±0.02 


66.5 ±1.0 


101.6 ±1.5 


IC131 


-0.06 ±0.03 


71.0 ±2.2 


107.2 ±2.4 


NGC588 


0.00 ±0.05 


69.8 ±2.6 


103.1 ±2.9 


IC142 


0.09 ±0.04 


77.8 ±4.6 


96.9 ±3.6 


B691 


0.27 ±0.05 


51.1±2.6 


93.0 ±3.4 


NGC592 


0.13 ±0.04 


64.5 ±2.6 


97.9 ±2.8 


M33 : 








R< 7.5 kpc 


0.72 ±0.04 


51.4 ±4.2 


63.2 ±5.3 



Using the standard method, iBuczilowski 



determined the scale length of the 6.3 cm ther- 
mal and nonthermal emission as 1.8 ± 0.2 kpc and 
4.2 ± 0.3 kpc, respectively. These scale lengths are 
smaller than those obtained here (even smaller 
than those obtained from the standard method), 
although the ratio of the nonthermal to thermal 
scale lengths is the same (—2). Due to the low 
signal-to-noise r atio of the o l d 6.3 cm Effelsberg 
receiver used by iBuczilowskil (|l98 8') much of the 
diffuse emission in the outer parts of the galaxy 
had been missed, leading to steeper radial profiles 
and smaller scale lengths. For the same reason also 
the sc ale lengths derived by iBerkhuiisen fc Klein 
( 1985f ) are too small. They obtained the distribution 
of the thermal emissi on at 6.2 cm from a catalogue 
of HII regions in Ha ( Boulesteix et al.iri974.) where 
the diffuse emission was not completely included. 

In case of equipartition between the magnetic 
field and CRs, the scale length of the CR electrons is 



given by /cr = ln(-i+an) /'2 and that of the magnetic 
field by = ^l^r (e.g. iKlein et al.lll982l ). Taking 
Q!„ ~ 1, we obtain Icr — 12 kpc and Ib — 24 kpc. 



For NGC6946. IWalsh et al.l (|2002f ) found a non- 
thermal scale length of In ~ 4 kpc which gives a 
smaller CR scale length of Icr — 8 kpc. Although 
NGC6946 is a Scd-type galaxy like M33, it is a star- 
burst system and it hosts stronger star formation 
in its central region than M33. This may cause a 
steeper radial profile of the total intensity and con- 
sequently a smaller nonthermal (and cosmic ray) 
scale length than in M33. 



9.3. North-south asymmetry 

From Figs. El and [6] it seems that the thermal emis- 
sion is stronger in the southern than in the northern 
half of M33. To investigate this north-south (N- 
S) asymmetry, we obtain the integrated flux den- 
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Table 5. Exponential scale lengths of the thermal 
and nonthermal components of the radio contiuum 
emission from M33. 



A 


Ith (kpc) 


l„ (kpc) 


new method 


20 cm 


2.4 ±0.2 


5.8 ±0.5 


3.6 cm 


2.6 ±0.2 


6.1±0.7 


standard method 


20 cm 


3.5 ±0.5 


6.2 ±0.7 


3.6 cm 


3.7 ±0.5 


8.9 ±0.9 



sity of the thermal, nonthermal, and total emission 
in each half separately. Table 6 shows the results 
from both the new and the standard method. The 
thermal emission from the new method is slightly 
stronger in the southern than in the northern half 
of the galaxy at both wavelengths. This may be 
the reason for the higher Faraday depolarization 
found in the south ern hal f than in the northern half 
2007b[) . In Paper III, we will dis- 
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cuss whether the asymmetry in the Faraday depo- 
larization can be caused by this N-S asymmetry in 
the thermal emission distribution. Note that, from 
the standard method, the thermal emission in the 
southern half is weaker. 

The nonthermal emission at 20 cm from the new 
method in the southern half is also stronger than in 
the northern half (in contrast to the 3.6 cm nonther- 
mal emission). This is expected as there are more 
supernova remnants in the southern half, ~ 56% of 
the total number, especially in the main southern 
arm IS (see Fig. |5]), which produce stronger non- 
thermal emission at longer wavelengths. The exis- 
tence of the strong optically thick HII regions (e.g. 
IC133) in the northern half of the galaxy also causes 
the different N-S ratios of the nonthermal emission 
at 3.6 cm. 

9.4. Uncertainties 

How do the assumptions of the new method of sep- 
arating thermal and nonthermal emission influence 
the results? We take the thermal fraction (e.g. at 
3.6 cm) as the final result and define U as the uncer- 
tainty in Ffii when one of the assumptions changes. 



U : 



th 



Fth 



(10) 



Fth and are the thermal fractions before and 
after changing an assumption, respectively. 

One of the assumptions is the choice of effec- 
tive extinction factor fd = 0.33 (non-uniform ioniza- 
tion). For a homogeneous distribution of dust and 



Te= lOOOOK 
Te= 8000K 
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Thermal fraction 



Fig. 14. Distribution (histogram) of the thermal 
fraction at 3.6 cm in M33. The number of bins used 
for this histogram is 50. 



ionized gas, fd= 0.5, the average uncertainty in the 
thermal fraction at 3.6 cm is only U — 2% with a 
standard deviation of 1%. 

Changing Tg from the adopted 10000 K to 
8000 K, the most probable uncertainty across the 
galaxy isU — 17% with a standard deviation of 3%. 
Fig'im shows that the thermal fraction decreases 
with decreasing electron temperature. Larger differ- 
ences in the thermal fraction are found for regions 
with higher thermal emission like HII regions, and 
the influence of changing Tg is relatively small in 
diffuse regions with small Fth- Thus, we find that 
the uncertainty resulting from the new method is 
mainly determined by uncertainty of the electron 
temperature. 

Another question is which method, the new 
method with its electron temperature assumption 
or the standard method with its constant nonther- 
mal spectral index assumption, faces larger uncer- 
tainty in the thermal fraction. For this comparison, 
we first obtain the nonthermal spectral indices for 
different electron temperatures. For =10000 K 
the most probable nonthermal spectral index is 0.9, 
while it is 0.8 for =8000 K. Then, from the stan- 
dard method, we calculate the thermal fractions as- 
suming a„ = 0.9 and 0.8, respectively, which leads 
to an uncertainty of U — 40%, much larger than 
that determined from our method (17%). We con- 
clude that the thermal fraction is more sensitive to 
variations of the nonthermal spectral index in the 
standard method than to variations of the electron 
temperature in the new method. 

Table 7 shows the thermal fractions obtained 
from the thermal and total flux densities of M33 
(integrated for R< 7.5 kpc) for different combina- 
tions of the assumptions of Tg and fd- 
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Table 6. North-South ratios of the integrated flux densities of the thermal, nonthermal, and total radio 
emission. 



A (cm) 








new method 


3.6 


0.89 ±0.05 


1.15 ±0.08 


1.10 ±0.08 


20 


0.86 ±0.05 


0.75 ±0.07 


0.82 ±0.07 


standard method 


3.6 


1.20 ±0.07 


0.67 ±0.07 


1.10 ±0.08 


20 


1.14±0.06 


0.77 ±0.08 


0.82 ±0.07 



Table 7. Average thermal fraction of M33 at 3.6 cm 
for different assumptions. The uncertainty U of 
each assumption is calculated with respect to the 
case Te = 10000 K and fd = 0.33. 
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10. Summary and conclusions 

We have developed a new method to separate the 
thermal and nonthermal radio emission from a 
galaxy. We used the highly resolved and sensitive 
Spitzer 70 and 160 iim. data of M33 to correct the 
Ha map of iHoopes fc Walterbod (|2000l) for extinc- 
tion. From this map, we calculated the thermal 
(free-free) emission at 3.6 and 20 cm and obtained 
maps of the nonthermal emission as well as a map 
of the nonthermal spectral index in M33. The dis- 
tribution of the nonthermal spectral index greatly 
helps to understand the origin and propagation of 
cosmic ray electrons in a galaxy. In brief, the results 
and conclusions are as follows: 



The distribution of the dust extinction is sim- 
ilar to that of the 160 /xm emission. The mean 
extinction in rings in the galactic plane exhibits 
a shallow radial gradiant. 

With a nonthermal fraction of about 30%-60% 
at 3.6 cm, the spiral arms and star- forming re- 
gions have a considerable contribution to the 
nonthermal emission. This contribution is negli- 
gible in the nonthermal maps obtained from the 
standard separation method. The radial profiles 
of the surface brightnesses and the wavelet spec- 
tra show that the distribution of the nonthermal 
emission from the standard method is smoother 
than that derived from the new method. This 



is caused by the assumption of a constant non- 
thermal spectral index in the standard method. 

• The nonthermal emission from the new method 
is still more smoothly distributed than the ther- 
mal emission. The exponential scale lengths of 
the nonthermal emission are more than twice as 
large as those of the thermal emission. 

• The standard method over-estimates the ther- 
mal fraction, especially at the position of 
giant HII regions. For galactocentric radius 
R<7.5kpc, the thermal fractions at 3.6 cm are 
51 ± 4% and 63 ±5 % from the new and standard 
methods, respectively. 

• For the first time, we derived a map of the non- 
thermal spectral index. In the star-forming re- 
gions, the nonthermal spectrum is relatively flat 
with an average value of of 0.6 ±0.1, the typ- 
ical spectral index of supernova remnants, but 
ttn increases to 1.2 ± 0.2 in the interarm regions 
and outer parts of the galaxy. This shows that 
the relativistic electrons lose energy when dif- 
fusing from their origin in star-forming regions 
towards interarm regions and the outer parts 
of the galaxy. The mean spectral index of the 
nonthermal emission becomes equal to that of 
the total emission at R~4.5kpc. This indicates 
that the total radio emission is mostly nonther- 
mal at R>4.5kpc in M33, where the spectral 
index is dominated by synchrotron and inverse 
Compton loss processes. 

• The north-south asymmetry in the distribution 
of the thermal emission obtained from the new 
method is opposite to that from the standard 
method. The thermal fractions are 47 ± 5% and 
53 ± 5% in the northern and southern half, re- 
spectively. 

• Generally, the integrated results from the two 
methods match with each other within the er- 
rors, indicating that the assumption of a con- 
stant nonthermal spectral index is a proper ap- 
proximation for 'global' studies. 
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